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CN ' Abstract 

Q>^ ■ A microscopic model for describing a superconducting mesoscopic weak link 

OO . 

^D . is presented. We consider a model geometry consisting of a narrow channel 

!>■ 

O 



coupled to wider superconducting electrodes which act as reservoirs fixing 
the asymptotic values of the complex order parameter. For this model, the 
H ■ Bogoliubov-de Gennes equations are discretized and solved self-consistently 
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using a non-equilibrium Green functions formalism. The transport properties 



a 
o . 

CJ . and the electronic excitation spectra of this system are studied for the different 



regimes that can be reached by varying parameters like coherence length, 
constriction length, normal transmission coefficient and temperature. We 
study in detail the transition from the point contact limit to the infinite 
channel length case, analyzing the maximum Josephson current that can be 
sustained by the weak link as a function of its transmission coefficient and 
length. It is also shown that for a constriction size ranging from zero to 
several times the coherence length, most of the current is carried, inside the 
constriction region, by bound states within the superconducting energy gap. 
These states correspond to Cooper pairs with binding energies smaller than 
the superconducting gap and which are spatially extended along the channel 
region, decaying exponentially inside the reservoirs. The importance of the 
self-consistent determination of the order parameter along the weak link is 



illustrated by analyzing different profiles obtained for channel lengths of the 
order of the coherence length. For temperatures not very close to Tc, our 
microscopic calculation predicts the appearance of features which cannot be 
obtained from Ginzburg-Landau theory. 

PACS numbers: 74.50.+r, 85.25.Cp, 73.20.Dx 

I. INTRODUCTION 

The interest in the physics of submicron superconducting devices [|I| stems to a great 
extent from the special combination, taking place in this kind of systems, of normal elec- 
tron phase coherence in the nanometer scale together with quantum macroscopic effects 
associated with the superconducting state. 

Some of the superconducting mesoscopic devices which are recently receiving more at- 
tention include normal metal-superconductor (NS) junctions and S-Semiconductor (Sm) 
junctions where Andreev reflection processes play a crucial role PH^, and superconducting 
mesoscopic point contacts which can exhibit quantization of the Josephson current . 

A theoretical description of the transport properties in all these submicron devices has 
to deal with spatial inhomogeneities due to the presence of geometrical boundaries and 
interfaces between different materials; in a superconducting mesoscopic device this could 
lead to strong spatial variations of the superconducting order parameter. 

In normal mesoscopic devices the problem of computing their transport properties has 
been historically addressed by different approaches. One is based on the scattering picture 
first proposed by Landauer and generalized to more complicated situations (multi-channel 
and multi-lead cases) by different authors |§|. Another approach, specially suited to deal 
with local inhomogeneities, relies on the use of a localized representation for the electronic 
states of the sample, whose transport properties can then be efficiently calculated in terms 
of Green functions PI. 



The scattering approach has been extended to superconductors by Blonder et al. |T^ 
and more recently applied to a mesoscopic NS junction by Beenakker [Q], who has obtained 
a multichannel generalization of the Blonder et al. result, and to different superconducting 
devices by other authors [0-0. Despite its many advantages, the scattering approach cannot 
describe properly situations in which the spatial variations of the superconducting order 
parameter on a scale smaller or comparable to the superconducting coherence length (^o) 
may be important. For instance, this can be the case in a NS junction when there is an 
appreciable induction of superconductivity in the normal electrode by a proximity effect 
||1 f II . Another example is that of a superconducting weak link with a length comparable to 
^0, in which case, as shown in the present paper, the self-consistent determination of the 
order parameter profile becomes unavoidable. 

In a recent letter |T2[ , we have presented a theoretical approach based on a local descrip- 
tion of the superconducting mesoscopic system, in which a non-equilibrium Green functions 
formalism is used. This method provides an efficient way of calculating the transport prop- 
erties and the spatial variations of the self-consistent order parameter through the whole 
system. The aim of the present paper is to discuss in further detail our approach, taking 
as a test system a mesoscopic superconducting weak link (SWL) at zero voltage. We shall 
consider a model geometry in which two wide superconducting electrodes are coupled by a 
narrow channel. As discussed by Likharev [l^ in his review on weak links, there are two 
main reasons for the basic interest on this kind of systems: first, weak links of reduced 
dimensions exhibit Josephson effect in non-tunnel conditions and are specially suitable for 
a variety of applications. On the other hand, for weak links of increasing length (larger 
separation between electrodes), the transition from the Josephson effect to bulk transport 
in superconductors can be studied. This last question is addressed in detail in the present 
paper. 

Traditionally, the transport properties of SWL have been analyzed with the help of 
phenomenological Ginzburg-Landau (GL) theory [|14|. However, GL theory is only valid 
for a limited temperature range (T -^ T^) and, furthermore, the implicit hypothesis of 



slow spatial variation of the order parameter in GL theory as derived from the microscopic 
theory |]15[, breaks down for a general mesoscopic geometry. In these situations a complete 



self-consistent solution of the Bogoliubov-de Gennes equations would be necessary |T^ . As 
we show in this work, features in the order parameter profile associated with length scales 
smaller than ^o, like the Fermi wavelength and some geometrical length scales, emerge in a 
natural way from the microscopic calculation. In addition, another result arising from the 
microscopic calculation, which could never be obtained from GL theory, is the existence of 
bound states inside the superconducting gap playing a fundamental role in the transport 
through the weak link. 

The plan of the present paper is the following: in section II we introduce our discretized 
(tight-binding) model for a general constriction geometry, and discuss the conditions for the 
fulfillment of current conservation within our model. In section III the non-equilibrium Green 
functions formalism in a superconducting broken-sjTiimetry representation is presented, giv- 
ing the expression for the relevant quantities within this formalism and details about its 
self-consistent determination. Section IV is devoted to an analysis of the transition from 
the point contact regime to infinite ID-superconductivity, discussing the maximum current 
that can flow as a function of the weak link transmission coefficient. The analysis of the 
spectral densities reveals the presence of a slowly increasing number of bound states inside 
the superconducting gap. It is shown that these states give the main contribution to the 
supercurrent in the constriction region and decay exponentially inside the 3D electrodes. 

In section V the importance of the self-consistent determination of the order parameter is 
illustrated by discussing the different type of profiles that can be obtained for a mesoscopic 
SWL. We emphasize those new features that cannot appear in a GL calculation. Finally, 
in section VI the effects of temperature both in the current-phase relationship and self- 
consistent profiles are analyzed. We find that, except very near Tc, there is no simple scaling 
of the system properties with temperature. 

The paper is closed with some concluding remarks. 



II. DISCRETIZED MODEL FOR A MESOSCOPIC WEAK LINK 

We consider a model weak link like the one depicted in Fig. 1. The constriction width, 
W, is assumed to be smaller than the penetration length and comparable to the Fermi 
wavelength, Xp, leading to a small number of conducting channels. The constriction length, 
Lc, can be varied from the point contact regime, Lc/C,o -^ 0, to the opposite case, Lc/C,o -^ cxd, 
recovering in this case the limit of a homogeneous quasi- ID superconductor. Phase coherence 
is assumed to be preserved along the whole system. 

The wider regions representing the left and right electrodes ensure the asymptotic con- 
vergence of the complex order parameter. A, to its bulk value. This model geometry would 



be similar to the ODSEE model proposed by Likharev |13 . 

Our aim will be the complete self-consistent solution of the microscopic Bogoliubov-de 
Gennes equations ||T^ for this system. For this purpose we find convenient to formulate 
these equations using a site representation for the electronic states. This representation can 
be viewed either as a tight-binding description of the electronic states or as a discretization 
of the Bogoliubov-de Gennes equations (a simple scaling of the parameters would account 
for the use of one or another description). For the electron-electron interaction we make the 



usual simplifying assumption of taking it as a contact attractive interaction II^, which in a 
site representation adopts the form of a negative-U Hubbard-like local term. 

Within these assumptions the mean-field model hamiltonian giving rise to the Bogoliubov 



de Gennes equations for our system can be written as fT^ : 



-n — 2^[^i f^jCia-Cicr + 2^ ^ij^ia^j^y 



JT^i.o- 



+ ^(A*c1^c1^ + A,qtqJ, (1) 

i 

where, for the zero voltage case, the chemical potential /x is a constant throughout the whole 
system and the sum over i and j in the second term is restricted to nearest neighbours only. 
The self-consistent conditions for the order parameter on each site are given by: 

A, = -U, < 4^4^ >, (2) 



where Ui is the attractive e-e interaction at site i. By choosing appropriately the f/j, tij 
and ej in the three regions (left electrode, constriction, right electrode) one can model dif- 
ferent situations: S-S'-S, S-N-S, S-Sni-S, etc. In this paper we shall concentrate in the first 
situation, assuming the same bulk value of the order parameter modulus, |A|, on the three 
separate regions and fixing the value ej = /i = (which implies e-h symmetry) in the whole 
system for simplicity. This criterion fixes the ratio f/j/tjj inside the three regions. The total 
phase drop along the whole system = 0^ — 0/j is imposed as a boundary condition, (p^ 
and (f)fi being the the bulk value of the order parameter phase on the left and right elec- 
trodes respectively. A relevant parameter is the superconducting coherence length in the 
constriction ^o(^)) which in our model can be estimated as ^o(^) = 2tc/7rA(T), where t^ is 
the hopping parameter in this region. 

The current between two neighbouring sites ij is given by: 

U^) = f E 0^- < 4At)c,At) > -tn < clit)cut) >) . (3) 

In the zero voltage case the supercurrent does not depend explicitly on time. It is worth 
noticing that current conservation is only fulfilled when the solution of the mean-field su- 



perconducting hamiltonian is fully self-consistent |T0|,|l2|,[T8|,|T9[ . Indeed, current conservation 
provides a stringent test of self-consistency. 

The proof of this statement is straightforward when using a site representation. Starting 
from the equation of motion for the electron density operator at site i {pi = X^o- Ha^^ia) we 
find: 

d < pi> i r. frl 

2ze 
= - E ^ij + l-(^i < c^TC^i > -A* < c^VIt >) (4) 

j '^ 

The last term appears due to the fact that H does not conserve the particle number. 
However, when the self-consistency condition (Eq. (2)) is fulfilled, this term vanishes, and 
the continuity equation, -^^^ + J2j hj = 0, is recovered for every site . 



III. SOLUTION IN TERMS OF NON-EQUILIBRIUM GREEN FUNCTIONS 

The averaged quantities appearing in Eqs. (2) and (3) can be expressed in terms of 
non-equilibrium Green functions [^ . For the description of the superconducting state it is 
useful to introduce spinor field operators ||2I|, which in a site representation are defined as: 



^i 






^J 



4t Cii 



(5) 



Then, the different correlation functions appearing in the Keldysh formalism adopt the 
standard causal form: 



C^fit^^t'p) = -t <T[i;,{t^)i;l{t'p)]> 



\(+' 



(6) 



where T is the chronological ordering operator along the closed time loop contour [^ . The 
labels a and /3 refer to the upper [a = +) and lower [a = — ) branches on this contour. The 
correlation functions Gti^ , which can be associated within this formalism with the electronic 



non-equilibrium distribution functions ||22|, are given by the (2x2) matrix: 
Gt-{t,t') = z 



( <c]^{t')c,^{t)> <9i(t')Qt(t)>^ 



(7) 



y< c]^{t')c\^(t) > <c,^{t')cl{t)> ) 

Eqs. (2) and (3) can then be written in terms of the Fourier transform matrix elements of 
Gtfit^t'): 
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For the zero voltage case the calculation of the different G^ (uj) is particularly simple 
because the following relation holds: 



Gt-{u;) = f{u;)[GUu^)-GUu;)] 



(10) 



where f{uj) is the Fermi distribution function, and G""'"^ are the advanced and retarded Green 
functions, which can be computed using recursive techniques [Q. (Note that this relation 
is the same as in a currentless state). 

The Green functions must be calculated self-consistently, according to Eq.(8). This is 
achieved starting from an initial guess for the order parameter profile and then using an 



iterative algorithm. As reported in [12|, the electrodes can be modelled in a simple way 
by Bethe lattices. This choice both facilitates the computation of the Green functions and 
ensures a fast spatial convergence of the order parameter to its bulk values. The Bethe 
lattice geometry is able to simulate in a simple way the geometrical dilution of the current 
taking place when passing from a quasi- ID to a 3D structure. We have fixed the Bethe 
lattice coordination number z + 1 = 4, which leads to a convergence of the order parameter 
within three or four layers. 

IV. TRANSITION FROM POINT CONTACT TO ID-FLOW 

It is interesting to analyze in detail the maximum d.c. Josephson current, Ic that can 
be sustained by a SWL with a single conducting channel as a function of its length and 
transmission coefficient, a. Within our tight-binding model a is a known function of the 
hopping parameters tij. For the results presented in this work a is varied by changing 
uniformly tij inside the constriction, while keeping them fixed inside the reservoirs. 

For any value of a, the point contact limit {Lc/C,o -^ 0) is well understood P,0,|^. In 



this case, the phase proffie can be well approximated by a step function and the current- 
phase relationship can be obtained analytically for a symmetric junction (details of this 
derivation within our formalism are given in Appendix A): 



where e(0) = ±|A(T)|a/1 — «sin^(0/2) denotes the position of bound states within the 
superconducting gap. These bound states give the main contribution to the supercurrent 



across the interface. As we shall see the existence of bound states inside the gap remains 
even for a SWL of length much larger than ,^o- It should be mentioned that an expression 
essentially equivalent to Eq. (11) was first derived by Kulik and Omel'yanchuk within a 
semiclassical approximation to the Bogoliubov - de Gennes equations |^ 

The maximum supercurrent in this point contact limit can be obtained as a function of 
a from Eq. (11). At zero temperature we find: 

/c,o(«) = ^(i-v/r^^) (12) 

Therefore the maximum possible current through a single quantum channel is eA/Ti 
whereas in the tunnel regime Icfl — eAa/2/i, which corresponds to the well known 



Ambegaokar-Baratoff value 7rA/2e-Riv [HI for the single channel case. The function Icfl{ci) 
is plotted in Fig. 2. Let us mention that the divergence in the derivative dlc^/da as a 
tends to 1 is a zero temperature feature, which is quickly smoothed at finite temperatures; 
at T = 0.2Tc, /c,o(tt = 1) is reduced by a factor of ~ 20%. This suggests that the experi- 
mental observation of the quantized critical current eA/h requires quite delicate conditions 
(i.e., perfect transmission, symmetric junction, very low temperatures, etc). 

At zero temperature, when Lc is increased keeping the value of a fixed, the maximum 
current can either decrease or increase from the L^ = value. In Fig. 3 we represent 
the behaviour of /(0) with increasing Lc for three values of the transmission a. Only the 
positive part of 1(0) is plotted. A general trend is the appearance of multivaluation for 
Lc/^o > 1; although the precise value for this threshold ratio is dependent on a, as can be 
observed in Fig 3. Note that the form of the /(0) curves itself is strongly dependent on 
a; therefore there is not a simple scaling between the different set of curves. For instance, 
the low transmission case (Fig. 3c) becomes markedly different from the other two when 
Lc/C,o increases (note the peculiar form of the 48-sites case in which there appears a cusp 
at = 27t; for larger lengths this cusp progressively bends down while /(0 = 27r) tends to 
zero). In the limit Lc/C,o — > oo, the behaviour would be that of a double tunnel junctions 



with a I{(j)) ~ sin(0/2) characteristic for the upper branch |]T2 



The maximuni current Ic is represented in Fig. 4 as a function of Lc/C,o for the three 
cases shown is Fig. 3 . It can be observed that when L^/^o ^ 1, and for high transmission 
Ic decreases asymptotically to a limiting value, whereas for low transmission the trend is 
the opposite. On the other hand, when Lc/C,o ~ 1 interference effects can lead to a non- 
monotonic behaviour of Ic with length. See for instance the intermediate transmission case 
(a = .75) where a dip at Lc ~ ^o is found. 

The limiting case Lc/C,o -^ oo corresponds to a homogeneous flow in a one dimensional 
superconductor in which the coupling to the left and right electrodes acts as boundary 
conditions. The asymptotic value of the maximum current Ic,od depends on this coupling 
and is therefore a function of a. These values are represented by the triangles in Fig. 2, 
where they can be compared with the Lc = case. It is worth noticing that for perfect 
transmission Ic^oo coincides with the depairing current of a one dimensional superconductor, 
i.e. 2eA/7Th [^. The derivation of this result within our model is given in Appendix B. 
The self-consistent phase gradient along the linear chain is equal to 2A/%vf when this limit 
is reached. 

On the other hand, for low transmission, when the matching to the electrodes is poor, the 
self-consistent phase drop concentrates mainly on the contacts and, as mentioned above, the 
system becomes equivalent to a double tunnel junction. In these conditions the maximum 
current is controlled basically by the transmission through a single junction. For a symmetric 
weak link this single junction transmission ai, is related to the total transmission a by 
«! ~ 2^0 (this relation holds only for a ^ 1). Then, according to Eq. (11) Ic,oo{ci) ~ x'V^) 
which is qualitatively in agreement with the numerical results for small a in Fig. 2. 

A deeper insight on the transition from the point contact to the infinite ID case can be 
obtained by analyzing the evolution of the local quasi-particle spectral density, Pi{uj) and 
the associated current density ji{uj), given by: 

p.ioo) = -Im[GUuj)]n (13) 

TV 

Mu) = ^ {u,+,[Gt\A^)]u - t,+i,,[G+^i(^)]n) (14) 

in 



where site i is chosen to be located at the center of the constriction, where the effect of 
the 3D reservoirs is less pronounced. In the following discussion only energies u < are 
considered as Pi{uj) = pi{—uj) due to e-h symmetry. 

For Lc/^o ^ 1 the most relevant feature is the appearance of a bound state in the spectra 
inside the gap. The position and weight of this state is a function of the phase difference, 
0, and the transmission coefficient, a, as noted in the discussion of Eq. (11). In this regime 
this state gives the main contribution to the current in the constriction: ji{uj) is essentially 
a delta function at the bound state energy. In Fig. 5 the evolution of Pi{uj) and ji{uj) as 
Lc/C,o increases is represented for the case of perfect transmission. In all the cases plotted 
in this figure the current is the maximum one for the corresponding constriction length. 

Two main features are observable in these curves: first, with increasing L^ new bound 
states appear in the gap. The new states initially split from the continuum, moving to 
energies closer to ;U = 0. Eventually, for L^. ^ oo these states would fill the gap as shown in 
the uppermost curves of Fig. 5, which correspond to the uniform infinite ID superconductor 
carrying the critical current (see Appendix B). The number of bound states increases very 
slowly with Lc/Co- For instance, for Lc/C,o as large as 10, only three bound states are present. 

On the other hand, the current density of the continuous part of the spectrum {uj < —A), 
which at Lc/C,o — ^ is negligible, becomes more important with increasing length, and in 
the limit Lc/^o ~^ oo gives a contribution which has the opposite sign to the total current. 

To conclude this section, let us analyze briefiy the nature of the bound states. The study 
of their spatial distribution shows that they correspond to Cooper pairs that are extended 
along the constriction region and decay exponentially inside the reservoirs within a typical 
length ~ ,^o- In Fig. 6 we represent the local spectral weight of the bound states along the 
constriction. It corresponds to case (d) in Fig. 5, in which three bound states are present. 
As can be seen in Fig. 6, the bound state closer to cj = has always a nodeless form, while 
the number of oscillations increases for the bound states with a larger binding energy. This 
situation is reminiscent of the one found for a potential well; however, one should keep in 
mind that the bound states in the present case are due to variations of the superconducting 

n 



phase on a finite spatial region (the constriction region), disappearing when the total phase 
drop is zero. In junctions of the type S-S'-S or S-N-S with large variations in the pairing 
potential one could find bound states even without current p|, p!7| , p8| . We should comment 



that the spatial distribution of the current carried by a bound state follows closely that of its 
weight. Actually, we have verified that the ratio between current and weight for each bound 
state is practically a constant along the constriction region. This result further illustrates 
the fact that the bound states correspond to Cooper pairs with a well defined velocity inside 
the constriction. 

V. SELF-CONSISTENT ORDER PARAMETER PROFILES AT ZERO 

TEMPERATURE 

The self-consistent determination of the order parameter profiles is essential for a weak 
link of length Lc> ^o, whereas for Lc '^ C,o, the point contact limit, the detailed form of the 
self-consistent profiles becomes irrelevant. In this latter case the current-phase relationship 
is given to a high degree of accuracy by Eq. (11) @,0]. 



In this section we present results for L^. > C,o in order to illustrate the effects of self- 
consistency. Although some of the overall features appearing in the profiles can be expected 
on intuitive physical grounds, their detailed form reflects a complex interplay between the 
different model parameters. 

In Fig. 7 we represent the phase and modulus profiles for fixed L^ and 0, and three 
different values of the transmission coefficient. All three results correspond to the upper 
branch of the /(0) characteristic. 

The common general features of the profiles along the upper branch are displayed in 
these figures. They consist on a constant phase gradient along the constriction together 
with localized drops at the contacts with the reservoirs. One could draw an analogy between 
these localized phase drops and the voltage drops at the contacts with the leads in a normal 
mesoscopic sample (Sharvin resistance) [§,^ . On the other hand, the modulus is on average 
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constant along the constriction. 

Superimposed to this general structure, oscillations with a spatial period Ai?/2 can be 
observed [0 . These oscillations have a maximum amplitude in the proximity of the contacts 
and decay when moving inside the constriction with a typical decaying length ~ ,^o- This 
behaviour is most clearly seen in case (c) of Fig. 7, which corresponds to Lc/^o ^ 1 and low 
transmission, leading to oscillations which are more concentrated near the contacts, having 
a larger amplitude due to the presence of larger "barriers" . This interference phenomenon 
is a consequence of the phase coherence of normal electrons along the mesoscopic system. 
As it is shown in the next section, there is a gradual disappearance of the oscillations with 
increasing temperature. Let us remark that a similar oscillation pattern can be found for 
the electrostatic potential along a normal mesoscopic constriction p9|-pl[l . 



When I{(j)) is a multivalued function, the solutions in the lower branch exhibit a very 
different character: they correspond to the solitonic-like solutions predicted by GL theory 
p!3| , p^j3^ . In Fig. 8 we show the phase and modulus profiles for this second type of solution, 
corresponding to three different values of L^ with fixed a and 0. We have chosen cf) r^ -k where 
the specific features of these solutions are more pronounced. It can be observed that the 
gross features are the ones predicted by GL theory, namely an abrupt phase drop together 
with a significant depression of the modulus which nearly goes to zero at the center of the 
constriction. Note that this profile leads to a very low value of the current. In addition 
to this general shape, an oscillation pattern of period Xp/'^ is again found. An interesting 
feature is the existence of a well defined "core" region of length ~ ^q where the phase drop 
takes place and the modulus is nearly zero (see Fig. 8). This core remains practically 
unchanged with increasing L^.. This is a feature that cannot be predicted by GL theory, 
the reason being, as we shall see in the next section, that the core size vanishes when T 
approaches Tc. 

The question about the stability of these type of solution deserves some attention. As 



discussed by Langer and Ambegaokar [^, these solutions correspond for a homogeneous 



infinite system to saddle points of the free energy. However, for a finite weak link these 
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solutions are presumably unstable |T3|. As a matter of fact, we have verified that these 
solutions cannot be reached using a simple iterative algorithm no matter how close from 
the actual solution one starts. Instead, the use of a more sophisticated algorithm (Broyden 
type [Q) capable of obtaining any sort of extrema, leads to the solitonic solution without 
difficulties. 

VI. EFFECTS OF FINITE TEMPERATURE 

In this section we discuss the effects of a finite temperature on the /((/)) characteristic 
and on the self-consistent profiles for our model weak link. Although L^/^oiT) is the main 
parameter controlling the different regimes (point contact, infinite homogeneous supercon- 
ductor) there is not a simple scaling of the system properties as a function of this single 
parameter. 

This becomes apparent when studying the evolution of the different properties with 
temperature. Fig. 9 shows the current-phase relationship for the same set of parameters 
as in Fig. 2b and three increasing values of temperature: 0.3, 0.5 and 0.95 T^. As one 
can observe, the phase interval where /(0) is multivalued gradually disappears while /(0) 
tends to the familiar sin(f) dependence as T — » Tc. This behaviour can be qualitatively 
understood noticing that for a given length and T sufficiently high (T -^ T^) , the condition 
Lc/C,o{T) ^ 1 is fulfilled, and the system should behave like a point contact in which case 
Eq. (11) approximately holds. Note, however, in Fig. 9c, that a significant deviation from 
the maximum current value predicted by Eq. (11), i.e Ic = ^Jr k t ^^ observed as soon 
Lc is a small fraction of C,o{T). This is due to the fact that the self-consistent phase profile, 
even for small Lc/^o{T), deviates from a step function form. 

The possible scaling of the maximum supercurrent with Lc/^o{T) can be analyzed with 
the help of Fig. 10. For these values of the transmission coefficient (a ~ 0.75), a universal 
behaviour of the type ~ exp~'^-^'=/'^o'^"^^ is only observed for temperatures larger than 0.5Tc. 
For lower transmission this departure from a universal behaviour is even more pronounced, 
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the maximum supercurrent at low temperatures is in this case an increasing function of 
Lc/^o(T) (see bottom curve in Fig. 4) while the exp^^^"^^"^'^^ behaviour is only recovered 
for T -^ Tc. 

Let us briefly comment the effects of finite temperature on the self-consistent order 
parameter profiles. Either in the upper and lower branch solutions there is a gradual disap- 
pearance of the Xf/'2 oscillations (see Fig. 11). In the upper branch solutions there are no 
notable changes in the phase profiles, besides the smoothing of the oscillations pattern, the 
expected Josephson current decrease being due to the global lowering of A(T) in the whole 
system. The evolution of the solitonic profiles with temperature is more unusual. This is 
illustrated in Fig. 11, where it can be observed how the solitonic "core", which at zero 
temperature has a size ^o(O), shrinks and eventually disappears as temperature increases. 
This evolution is particularly clear in the phase profiles of Fig. 11. 

VII. CONCLUDING REMARKS 

We have presented a self-consistent solution of the microscopic Bogoliubov-de Gennes 
equations for a mesoscopic SWL. As illustrated in the present work, our method allows 
us to analyze situations where spatial variations of the superconducting order parameter 
over length scales smaller or comparable to the coherence length are important, and its 
self-consistent determination becomes necessary. As discussed in this work, the usual GL 
approach may not be valid for a mesoscopic sample at low temperatures. Our results reveal 
that for low temperatures and due to coherence effects, there is no universal behaviour of 
the SWL properties as a function of Lc/^o{T). These coherence effects are also reflected in 
the self-consistent order parameter profiles, which exhibit features that cannot be predicted 
by GL theory. 

We have analyzed in detail the transition from the point contact limit to the infinite 
ID superconductor characterized by a continuous spectrum. The appearance of bound 
states seems to be a general feature of a weak link of mesoscopic size due to the spatial 
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inhomogeneity of the pairing potential in the constriction region. Even in the absence of a 
well defined pairing potential well, spatial variations of the superconducting phase associated 
with a supercurrent can lead to the formation of bound states. 

The application of the present approach to the description of submicron superconducting 
devices of current experimental interest |l| , which may include S-N or S-Sm interfaces and 
where the sample geometry plays an important role, is under progress in our laboratory. 
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APPENDIX A: 

In this appendix we give a short account on the derivation of the supercurrent-phase 
relationship for a SWL in the point contact limit. Starting from Eq. (9), the current 
evaluated at the interface between the electrodes can be written as: 

Op /"oo 

I=-f^j_^ dujf{u) Tr [TUG^L - Gk) - T^nUGln - G^^)],, (Al) 

where L and R indicate the left and right electrode and Tr denotes trace over orbitals at the 
interface. The hopping elements in the superconducting broken-symmetry representation 
are given by 2 x 2 matrix Tjj = r^tij, where T3 is a Pauli matrix and i,j denote any pair of 
orbitals. 

The Green functions appearing in this expression can be obtained from the ones corre- 
sponding to the uncoupled electrodes (Tlr = 0) g^'*^, g^*^, using Dyson equations: 

Gl^ = glT.Rg'iPl^ (A2) 

where D^/j = [I — T^/jg^T/jigJ]" with similar expressions for the retarded quantities. 
Substituting in Eq. (Al) and performing some elementary algebra, we obtain: 



I=^j_^ du^fico) Tr [iB%,, - D2^) - (D^^ - D^^)],, (A3) 

This equation is valid for a general junction geometry with any number of conducting 
channels. For the simple case of a symmetric junction with a single channel, an analytical 
expression can be derived for J(0) if a step-like phase-profile is assumed: 



T 4e« 2 . , /""^ > ./ N 91,2iM9r,i2{^) 9l,2i{^)9r,12{^ 
^ = -r'^LR sm0/ dujf{uj) 



(A4) 



where the tilde indicates that the phase factor has been removed, i.e. gL,2i = exp"**^^ gL2i 
and Qr^u = Q^V^'^^ 9r,i2- The integration in Eq. (A4) can be performed analytically as 
a contour integration by realizing that the main contribution (up to corrections of the 
order A/ep) is given by the zeroes within the superconducting gap of det{D1^{u!)) , which 
correspond to the bound states commented in section IV. The contribution from these poles 
yields Eq. (11) straightforwardly. 

APPENDIX B: 

A simple derivation of the depairing current for a ID superconductor within our tight- 
binding model is given below. 

The hamiltonian in this case is the one given in Eq. (1) with tij = t{Sij+i + Sij-i), 
ej = and Aj = Ag exp*''-'. This last relation holds for a uniform superconducting flow, with 
Cooper pairs of net momentum q |T^ . For an infinite system, the Green functions adopt a 
simpler form in the ^-representation: 

[Gm(^)]ii = 7 ,7 tu'"'^' ^ A2 (Bl) 

{uj - ek){u + e_k+q) - A^ 

[GmM]i2 = -7 Tr4^ r^r^ ^^2) 



where ek = 2tcosk — fi. The poles of these functions at Ej^^ = Sek,q ± J^l^q + A^, where 
6ek,q = (cfc — e_fc+g)/2 and ek,q = (e^ + e-.k+q)/'^, give the excitation spectrum of a current 
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carrying ID superconductor. The critical momentum q^ which corresponds to the condition 
that the gap in the excitation spectrum goes to zero, is given, for sufficiently small qc 
{qcSinkp <^ 1), by gc = Ao/t sin kp. 

It can be shown that the self-consistent order parameter A^ for q < qc is in this limit 
equal to Aq up to corrections of order g^. Actually, these corrections are positive, leading 
to a small increase in A^ with increasing current. 

The local spectral density p{uj) and the local current density j{uj) discussed in section 
IV can be directly calculated from Gk,q{uj). In particular, the total current is given by: 

I=—t dksink duf{uo)Im GlAuo) (B3) 

At zero temperature, the integral over u can be easily performed giving: 

et f^ ,, ek^qsink 2e 



^ , dk r =^^tsinkF q (B4) 

TihJ-n . el +A2 Tih ^ 

y k,q r ^q 

where the last approximation holds in the limit qsinkp ^ 1. Thus, for q = qc the depairing 
current is simply -eA/h, in agreement with ref. ||27|| . 
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FIGURES 

FIG. 1. Schematic representation of our discretized model weak link. 

FIG. 2. Maximum supercurrent for a single channel weak link at zero temperature in the point 
contact limit {Lc = 0) and in the asymptotic infinite length limit, as a function of the transmission 
coefficient. 

FIG. 3. Supercurrent-phase characteristics at zero temperature for different constriction lengths 
and three pair of values for the transmission and the coherence length: (a) a = 1 and ^q = 22.05; 
(b) a = .75 and ,^o = 12.73; (c) a = .28 and ^q = 6.37. The coherence length and the Lc values 
given in the figure are measured in units of the site spacing. 

FIG. 4. Maximum supercurrent at zero temperature for the same three cases as in Fig. 3 as a 
function of Lc/^o- 

FIG. 5. Local spectral and supercurrent density at the center of the constriction for the a = 1 
case and different values of Lc/Co'- (a) Lc/^o = 2.18, (b) Lc/$,o = 5.80, (c) Lc/$,o = 9.43, (d) 
Lc/Co = 13.06 and (e) Lc/^o = oo. The curves are displaced upwards with respect to case (a) for 
clarity. All cases correspond to the maximum supercurrent. 

FIG. 6. Spectral weight along the constriction for the bound states appearing in case (d) of 
Fig. 5. The order from top to bottom corresponds to increasing binding energy. 

FIG. 7. Upper branch self-consistent order parameter profiles (phase and modulus) for a total 
phase drop (j) = 4.60 and fixed constriction length Lc = 88. The transmission and coherence length 
correspond to the three cases shown in Fig. 3. The curves (a) and (b) are displaced upwards with 
respect to case (c). 

FIG. 8. Lower branch self-consistent order parameter profiles (phase and modulus) for a total 
phase drop (j) = 3.20, a = .75, ^o = 12.73 and three values of the constriction length: (a) Lc = 48, 
(b) Lc = 64 and (c) Lc = 88. 
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FIG. 9. Supercurrent-phase characteristics for a = .75, different constriction lengths and three 
values of temperature: (a) T = 0.3Tc and ^o(r) = 12.78; (b) T = CST^ and ^o(^) = 13.37; (c) 
T = 0.95Tc and i^iT) = 39.79. 

FIG. 10. Maximum supercurrent for a = .75 as a function of Lc/S,oiT) for increasing values of 
temperature. 

FIG. 11. Lower branch self-consistent order parameter profiles (phase and modulus) for a total 
phase drop (p = 3.20, a = .75, Lc = 88 and three values of temperature: (a) T = 0, (b) T = 0.3Tc 
and (c) T = 0.5Tc. 



9S 



